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ABSTRACT 

We present precise H I 21 cm absorption line redshifts observed in multiple epochs to directly constrain 
the secular redshift drift z or the cosmic acceleration, Av/Afo. A comparison of literature analog spectra to 
contemporary digital spectra shows significant acceleration likely attributable to systematic instrumental errors. 
However, we obtain robust constraints using primarily Green Bank Telescope digital data. Ten objects spanning 
7 = 0.09-0.69 observed over 13.5 years show z= (-2.3 ± 0.8) x 10"^ yr"' or Av/Afo =-5.5 ± 2.2 m s"' yr"'. 
The best constraint from a single object, 3C286 at (z) = 0.692153275(85), is z = (1.6 ± 4.7) x lO"** yr"' or 
2.8 ± 8.4 m s"' yr"' . These measurements are three orders of magnitude larger than the theoretically 



Av/Af, 

expected acceleration at z = 0.5, z = 2 x 10~'' 



yr or Av/A?o =0.3 cm s yr , but they demonstrate the 



lack of peculiar acceleration in absorption line systems and the long-term frequency stability of modem radio 
telescopes. A comparison of UV metal absorption lines to the 21 cm line improves constraints on the cosmic 
variation of physical constants: /S.(a^ g p n) / a~ g p ^ = (-1.2 ± 1.4) x lO"*" in the redshift range z = 0.24-2.04. 
The linear evolution over the last 10.4 Gyr is (-0.2 ± 2.7) x lO"'*" yr"', consistent with no variation. The 
cosmic acceleration could be directly measured in 125 years using current telescopes or in ~5 years using a 
Square Kilometer Array, but systematic effects will arise at the 1 cm s"' yr"' level. 

Subject headings: cosmological parameters — Cosmology: observations — Cosmology: miscellaneous — 
dark energy — quasars: absorption lines 



1. INTRODUCTION 

The measurement of secular redshift drift can directly de- 
termine the cosmic acceleratiorQ and the history of expansion 
in a model-independent manner (iSandagdl 1 9621: iLoebll 1 9981) 
and has recently been explored in the contex t of dark energy 
and large aperture opti cal telescopes (e.g., iCorasaniti et alj 
I2007t iLiske et al.ll2008h . The magnitude of the redshift drift 
is minuscule, of order 5 cm s"' decade"' for a galaxy at red- 
shift z = 1, simila r to gravitational acce lerations in galaxies 
and clusters (e.g.. lAmendola et alj|2008l ). Measuring this ef- 
fect requires extreme redshift precision over time baselines of 
decades, and the signal may thus be overwhelmed by other 
systematic effects, both physical and observational. For ex- 
ample, secular redshift drift may be induced by proper ac- 
celeration of the observer. Such proper acceleration (but 
not the associated redshift drift) has been very precisely 
measured using pulsar timing in a Galactic reference frame 
jZakamska & Tremaind l2005b and in a cosmological con- 
text usi ng the proper motion of extragalactic radio sources 
(iTitov e t al. 201 1; Xu et al. 2012). Proper acceleration of the 
solar system barycenter about the Galactic center creates an 
extragalactic dipole proper motion secular aberration drift sig- 
nature on the sky that distinguishes it from a cosmological 
signal: objects appear to be moving toward the Galactic cen- 
ter. Astrometry and pulsar timing cannot, however, detect a 
change in the expansion rate of the universe, which is purely 
radial. 

The observed rate of change of redshift, Az/Afo (hereafter 
z), is the difference between a constant expansion rate at red- 
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^ We use "acceleration" in the physical sense, which includes deceleration. 



shift z, Ho(l +z), and its actual value H{z): 
^=Ho(l+z)-H(z), 

where 



(1) 



H{z) = //o ^ r!M,o ( 1 + z)3 + 1] A + ( 1 - r!M,o - Oa)( 1 + z)2 (2) 

for a matter- and cosmological constant-do minated uni verse, 
and Ato is the observer's time increment (lLoeblll998h . For 
a time-varying dark energy parameterized with scale factor- 
depen dent equation of state w(a) = Wo + Wa(l-a) (iLinded 
12003b where a= 1/(1+ z) is the scale factor, the Hubble ex- 
pansion becomes 



Hiz) = Ho\ flM.oil +z)3 + fl„o(l +z)3(l+'-o+H'=)e-3»V,z/(l+z) 



+ (l-r!M,o-a.,o)(l+z)2. (3) 



The observed acceleration is 

Av _ Az/Afo 

aTo'"' (1+z) ■ 



(4) 



The secular redshift drift and correspon ding acceleration are 
plotte d for several cosmologies in Fig.fTt'Ba lbi & Ouercellinil 
120071) explore less conventional cosmological models and 
theories of gravity. A direct measurement of z is completely 
model-independent and directly indicates the history of ex- 
pansion, whatever that may be. 

Here we examine the feasibility of measuring or constrain- 
ing secular redshift drift using H I 21 cm absorption line sys- 
tems. We demonstrate that with modern telescopes, observa- 
tions spanning roughly a decade can measure redshifts to parts 
per 100 million and constrain the redshift drift to down a few 
m s"' yr"'. Analog spectra show a systematic redshift offset 
from modern digital spectra, which limits the time baselines 
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Fig. 1 . — Theoretical acceleration vs. redshift. Top: Secular redshift drift (Eqn.[T). Bottom: Apparent acceleration (Eqn.|4). The solid lines show the cosmology 
Ha = 12 km s"' Mpc"' , SIa = 0.73, and Q.m = 0.27, the dash-dotted lines show this cosmology with a varying dark energy equation of state, w(a) = Wo +Wa(\-a) 
(Lindei 2003), with n'o = -1 and hv, = ±0.5 (lower and upper curves, res pectively), the long-dashed lines show a matter-only universe with Q.m = 0.27 (open), 
and the short-dashed lines show a closed matter-only universe (Qm = 1.0). IBalbi &"O uercellmi l200l examine i for a panoply of less conventional models. The 
horizontal dotted lines indicate stationary redshift and velocity (crossing the concordance cosmology line at z.~q = 2.48), and the vertical dotted lines show the 
matter-cosmological constant equivalence redshift (zma = 0.39) and the transition from a decelerating to an accelerating universe (zg=o = 0.76). 



available forz measurements and impacts previous constraints 
on the cosmic evolution of physical constants that employed 
analog 21 cm spectra. We revise the measurement of the evo- 
lution of a^gpine/mp over the redshift span z = 0.24-2.04. 

In the following treatment of z, we make no assumptions 
about cosmology; the only numbers needed are the speed of 
light, c = 299792.458 km s"', and the rest frequency of the H I 
21 cm spin-flip transition, i^oiHI) = 1420.405751768 MHz, 
although neither is truly needed to measure a redshift change 
since it is ultimately simply a shift in the observed frequency. 

2. OBSERVATIONS, DATA, AND DATA REDUCTION 

Analog data provide lo ng time baselines , in some cases 
more than 30 years (e.g. IWolfe et anil978h , but we find a 
systematic offset toward lower redshifts compared to single- 
telescope digital epochs (Figure |2]i. These systematic offsets 
appear in most pre-digital spectra but are not consistent be- 
tween absorption line systems, suggesting that they are instru- 
mental and uncorrectable. We thus excluded the high qual- 
ity and potentially useful observations from analog epochs. 
Moreover, most extant digital spectra lack the signal-to-noise 
needed for this study, and in order to minimize systematics we 
restrict all spectral epochs to those obtained using the Green 
Bank Telescop43 (GBT), except as noted below. 

We observed 13 H 1 21 cm absorption line systems in 2003- 

' The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



2004 using the GBT Spectrometer and repeated the obser- 
vations in 2012 (programs GBT 03C-009 and 12A-134, re- 
spectively). Five of the 13 sources include a third epoch ob- 
tained at the GBT (two) or elsewhere (three). We reduced 
the archival GBT Spectral Proce ssor spectra of 023 5-H164 and 
3C286 (program GBT 07 A-02 1 : IWolfe et al.l20TTl) in an iden- 
tical manner to all other GBT spectra, d escribed below. For 
early third epochs (ca. 2000), we used the lCarilU et all (l2000h 
European Very Long Baseline Interferometry Network spec- 
trum of B0218H-357, the Kanekar & Chengalur ( 2001) Gi- 
ant Metrewave Radio Te lescope time-av e raged spectrum of 
PKS 1127-145, and the iLane & Bnggsl (|200T1) Westerbork 
Synthesis Radio Telescope spectrum of P KS 1243-072 (all 
digital observations, extracted by DEXTER: iDemleitner et al.l 
(jlOOl)). We observed 5 additional absorption line systems in 
single epochs (TablefTli. 

The two main observational epochs described here span the 
transition from analog to digital television. Consequently, ab- 
sorption line systems that were observable through the analog 
radio frequency interference (RFI) comb in 2003, particularly 
those at z > 1, were no longer observable amid the digital RFI 
in 2012, and lines that were lost to RFI in 2003 at the upper 
end of the old analog television allocation, PKS1629H-120 and 
3C216, were detectable in 2012. These single-epoch observa- 
tions are listed in Table[T]and are considered in our treatment 
of fundamental physical constants but are not used in the anal- 
ysis of Z- 

GBT observations were made in 5-minute position- 
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TABLE 1 

Redshift and Acceleration Measurements 
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ZHI 
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PKS 1243-072 
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1 


3 


56136 


4343 
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122(107) 


0235+164 


0.523741603(72) 


4 


3 


55953 


2699 


-1.3(1.5)xlO"' 


-25.9(28.7) 


PKS 1629+120 


0.5317935(11) 
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1 


56136 








B3 1504+377 


0.67324197(79) 


1 


2 


55968 


2978 


2.3(2.0) xlO"' 


40.6(35.3) 


B0218+357 


0.6846808(13) 


1 


3 


56038 


4936 


2.1(2.6)xl0"' 


37.1(46.2) 


3C286 


0.692153275(85) 


1 


3 


55947 


2956 
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2.8(8.4) 


PKS 1830-211 


0.8848633(41) 
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2 


56039 


3048 


-1.2(1.0)xl0-* 


-192(163) 


1331+170 


1.7763904(66) 


1 


1 
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1157+014 


1.943670(12) 


1 


1 


53070 








0458-020 


2.0393767(21) 


2 


1 


53071 









Note. — For each absorption line system, zhi is the error-weighted mean redshift of the strongest com- 
ponent in the barycentric reference frame, Nurxi. is the number of components used to determine z, Atibs is 
the number of observations, to is the Modified Julian Date (MJD) of the most recent observation, A/o is the 
full time span of the observations, ; is the secular redshift drift, and Av/A(o is the acceleration (Eqn. |4). 
Parenthetical values are \a uncertainties. Objects in bold are those used in the i analysis and Figure[3l those 
in italics are poorly constrained due to broad blended lines (3C196 and PKS 1830-21 1) or low signal-to-noise 
(PKS 1243-072). We include objects with a single digital observation for posterity, and some are compared 
to UV absorption fines in Figure|4] 

switched (nodding) scans using two linear polarizations, 9- 
level sampling, 2-3 s records, a 1-1.5 s winking calibra- 
tion diode injection signal, and a 12.5 MHz bandwidth cen- 
tered on the H I line. Typical on-source integration times 
were 1-2 hours. We subtracted the corresponding off-source 
record from each individual on-source spectrum and normal- 
ized the difference spectrum by the sky spectrum. Each spec- 
tral record was manually examined and flagged for RFl then 
averaged in time and polarization. After polynomial baseline 
subtraction and Hanning smoothing, we typically achieved 
3-24 mJy rms noise in 1.1-1.6 km s"' channels at 1300- 
750 MHz {z = 0.09-0.89) in 2003-2004, and 6-30 mJy rms 
noise in 0.35-0.61 km s"' channels in 2012 over the same 
redshift range. We used GBTlDlfl for aU GBT data reduc- 
tion. 

Absorption lines were fit using Gaussian profiles in each 
epoch and were not constrained by other epochs. Residual 
spectra were usually indistinguishable from line-free noise. 
Typically, only high signal-to-noise narrow well-separated 
lines were used for the z measurement, and Table [T] lists 
the number of components used in each case. We measured 
the error-weighted mean redshift of the strongest Gaussian- 
fit component in each system to provide a reference redshift 
for posterity and to indicate the degree of precision attained 
(TablelB. 

Figure |2] shows the spectra and time series of 0235H-164, 
including a pre-digital epoch, as an example of the statisti- 
cal uncertainties and systematic variation of the absorption 
lines. Weak components in the wings are only needed for 
good fits when the signal-to-noise is high, so these compo- 
nents are not included in the z measurement. We subtract the 

GBTIDL ( http : / / gbtidl ■ nrao ■ edu/ ll is tfie data reduction pack- 
age produced by NRAO and written in the IDL language for the reduction of 
GBT data. 
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Fig. 2. — Example of redshift drift observations of an H I 21 cm absorption 
line system, 0235+164. Top: Spectra obtained at the GBT (upper three, all 
digital) and Arecibo (analog, bottom; Wolfe et al. ( 1978)) labeled by MJD, 
scaled to match the MJD 55953 blue line, and offset for clarity. Gaussian 
fits to well-defined lines are color-coded (others are black), matching points 
in the lower panel. As with many of the absorption fine systems listed in 
Table [T] the analog spectrum of 0235+164 is systematically offset from the 
digital GBT spectra and is thus omitted from the linear fit. Bottom: Redshift 
(upper) and redshift mean offset (lower) versus MJD. Dotted lines indicate 
the no-drift (mean) locus in both panels, and solid lines show the linear fits, 
omitting the MJD 43032 analog epoch. Dashed lines show the formal error 
on the slope from an en'or-weighted linear least-squares fit. 
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Fig. 3. — Acceleration vs. redshift of H I 21 cm absorption line systems. Horizontal lines indicate the error-weighted mean (solid) and its Icr error (dashed). 



The top plot shows z: the bottom plot shows acceleration Av/ Ato. The error- weighted means are (z) -- 
which are dominated by 0738-1-313 A, 0738-H313 B, PKS 1413-1-135, and 3C286. 



(-2.3 ± 0.8) X 10-"* yr-' and (v) =-5.5 ± 2.2 m s"' yr" 



error-weighted mean redshift from each component and then 
make an error-weighted Unear fit to all reliable components 
simultaneously in order to average out stochastic variation 
among the components and obtain a bulk redshift drift of the 
absorber line ensemble. The slope of the linear fit is z- We es- 
timate the uncertainty in z from a Monte Carlo process except 
for large A^unes x A^obs systems 0738H-313 A, PKS 1 127-145, 
0248-1-430, and 0235-1-164, which use a formal error-weighted 
least-squares fitting error. The acceleration, Av/Afo, is cal- 
culated via Equation |4] and included in Table[T] 

3. RESULTS AND ANALYSIS 

Figure [3] shows the secular redshift drift and acceleration 
of 10 absorption line systems. While the theoretical z curves 
shown in Figure [T| could be fit to these data and would sug- 
gest a roughly linear trend for this redshift range, the typical 
measured z value is about three orders of magnitude larger 
than the theoretical expectation and is consistent with zero, 
so we simply calculate the mean z across redshift to assess 
the precision of this method and the presence (or absence) of 
systematic offsets from the expected value of zero. 

We measure an error-weighted mean secular redshift drift 
of {£) = (-2.3 ± 0.8) X 10"^ yr"' and an acceleration of 
(v) = -5.5 ± 2.2 m s"' yr"', both consistent with zero. An 
error-weighted linear fit fixed to z = at z = has slope 
-3.2 ± 13.3 m s"' yr"' per unit redshift. These measurements 
agree with expectations and show no evidence of secular red- 
shift drift or other systematic effects that might impair efforts 
to obtain higher precision, either via longer time baselines, 
higher signal-to-noise, or a larger sample (Sec. 14. 11 1. 

Given the systematic analog-digital offsets in HI 21 cm 
lines described in Sec. |2] it stands to reason that compar- 
isons of literature analog H I lines to other absorption lines 
in order to measure physical constants at various redshifts 
are systematically affected. The redshift difference between 
rest-frame UV metal and H 1 21 cm absorption lines measures 



(or constrains) the cosmic evolution of gpnie/nip, where 
a = /Tic ~ 1/137 is the fine structure constant, gp is the 
proton gyromagnetic factor, an d nig and nin are the electron 
and proton mass, respectively (IWolfe et al.lll976^ . Since the 
ratio of the H I 21 cm line to the UV metal lines scales as 
cP- gp fJ,, where /i = nig /nip, the redshift difference between the 
lines can be related to fundamental constants as 



A(a^ gplj) _ zuv-zii 



(a^gpfi)c 



1+Z2I 



(5) 



p rovided the gas is veloci ty-coincident. 

iTzanavaris et al.l (l2007h employed a number of analog mea- 
surements of the 21 cm lines presented in Table [T] to examine 
the cosmic evolution of gp /i, so we revisit their analysis us- 
ing digital data. Figure|4]shows the normalized redshift differ- 
ence between H 1 2 1 cm lines and U V metal lines for the nine 
absorption line systems measured bv ITzanavaris et akl (l2007h 
and for the seven systems measured by us usi ng the GBT. The 
UV lin e ensemble redshifts are presented in ITzanavaris et alJ 
(l2007h . In two cases we use a 21 cm component that is not the 
dominant one for comparison to the UV lines: in 0235-1-164 
we use the H I line at z = 0.52384114(41), the strongest of 
the line triplet (Fig. |2l ), to match the component used by 
ITzanavaris etaP (12007^ . and in 0458-020 we use the redder 
component at z = 2.0395474(64) that ma tches the many UV 
lines presented in lTzanavaris et al] (120071) . 

Despite the smaller sample, we obtain a substantially re- 
duced error- weighted mean: the fractional change in gpfjL 
is (-1 .2 ± 1 .4) X 10"^ in the redshif t range z = 0^ 4-2.04. An 
identical treatment of the Tzanavari s et al.l ( l2007h data yields 
a fractional change of (-1-19.9 ± 3.0) x 10"^ (using their statis- 
tical error bars for consistency). The unweighted means are 
simil ar, (-1-2.9 ± 5.5) x lO"'' ( this work) versus (-1-5.3 ± 9.9) x 
10"^ ( ITzanavaris et al.ir2007l) . and our uncertainty is roughly 
half despite the slightly smaller sample. 
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Fig. 4. — The fractional redshift difference between UV metal and 21 cm H I absorption lines vs. redshift. This is equivalent to the fractional change in 
gp fi. Top: Reproduction of the lTzanavaris et al., 1,2007) plot using literature 21 cm absorption Hnes and their statistical eiTor bars. Bottom: the same plot using 
GBT-measured H I 21 cm lines. Horizontal lines indicate the null line (dotted), the mean (red, short dashed), and the error-weighted mean (blue, long-dashed). 
Error bars in the lower left of each panel indicate the uncertainties in each mean value. Despite the smaller sample, our mean offsets are significantly smaller and 
their uncertainties are smaller by a factor of two. Individual error bars are likewise substantially smaller due to higher signal-to-noise in most 2 1 cm observations. 



The linear evolution over the last 10.2 Gyfl assuming the 
laboratory value at z = 0, is (-0.2 ± 2.7) x 10"'^ yr"', consis- 
tent with no variation. A negative value indicates a positive 
slope in Fig. |4] (the value decre ases as time advances f rom 
the early universe to the present). iTzanavaris et al.l (l2007h ob- 
tained a limit one order of magnitude larger, (-0.6 ± 1.2) x 
10"'^ yr"', using their favored analysis. 



4. DISCUSSION 
4.1. z 



because they will be randomly distributed, showing both pos- 
itive and negative values, and will not be coherent across the 
sky or depend on redshift. The barycenter acceleration will 
appear in z measurements as a redshift-independent accelera- 
tion dipole on the sky. 

Our z measurement depends primarily on the four best- 
measured absorption line systems; omitting points with error 
bars greater than 20 m s~' yr"' has a negligible impact on the 



result, which suggests that only 



10 additional systems with 
would improve this 



The intrinsic stability of the H I 21 cm lines does not seem 
to be an issue at the level of a few m s"' yr"'. Most sys- 
tems show negligible epoch-to-epoch redshift variation, but a 
few systems such as PKS 1127-245 do show significant ap- 
parently stochastic variation, and the measured z may be af- 
fected by the time of measurement and the number of epochs. 
Stochastic line shifts of order 500 m s"' over time scales of 
about 1 year with no net change in the mean line velocity 
have been observed (Wolfe et al. 1982; Briggs 1983) as well 
as variation in line depth (e.g., lKane kar & Chengalur 2001), 
and the common interpretation is that a variable illuminating 
radio source structure picks out varying portions of the fore- 
ground absorbing clouds. This effect may become significant 
as one approaches the precision needed to detect the expected 
cosmological acceleration. 

Gravitational accelerations of order 1 cm s"' yr"' are likely 
in galactic and extragalactic settings, such as the accelera- 
tion of g alaxies in clusters and th e orbital accelerations within 
galaxies (lAmendola et alj|2008l) . For example, the Solar Sys- 
tem barycenter accelerates toward the Galactic Center at ^^0.7 
cm s"' yr"', assuming a circular orbit. Extragalactic gravi- 
tational accelerations are distinguishable from the cosmic z 

^ From present to z = 2.04, assuming Ho = 72 km s"' Mpc"' . = 0.73, 
and Qm = 0.27. 



uncertainties of order a few m s ' yr ' 
measurement substantially, perhaps below 1 m s"' yr"'. 

To estimate the time required to measure the canonical cos- 
mological z at z = 1 (£ = 2 X 10"" yr"'), we assume a scenario 
whereby an ensemble of 30 lines is observed once per year 
with uncertainty a~ = 5 x 10~^ per line per observation. We 
claim a detection when the uncertainty of a linear fit of red- 
shift versus time recovers the slope with at least 3<t signifi- 
cance. For the GBT, a direct detection of z will require about 
300 years. For a Square Kilometer Array (SKA) with ^^130 
times the collecting area, CT; ~ 4 x 10"'" and this measurement 
can be made in about 12 years, which is competitive with a 
42 m o ptical telescope ob serving Lyman a forest lines over 
^20 yr (iLiske et al.ll200 8f). If the number of lines increases to 
300, an expected dividend of upcomin g SKA prototype sur- 
veys such as those p lanned for ASKAP (iJohnston et al.ll2008t 
iDarling et al.l 1201 lb , then the GBT measurement would re- 
quire ^^125 years, and the SKA measurement ^-^5 years. 

4.2. Physical Constants 

H I 21 cm lines are no longer the Umiting factor in con- 
straining the cosmic evolution of gpji: UV metal line red- 
shift uncertainties are 1-2 dex larger than the 21 cm line un- 
certainties at z < 1 and are comparable at z > 1 . Systematic 
effects, such as the coincidence of the 21 cm- and UV line- 
producing regions do not yet appear to be significant. Sys- 



6 



Darling 



tematic offsets between analog and digital 21 cm redshifts 
suggest that other astronomical tests of the cosmic variation 
of physical constants employing analog 21 cm data should be 
re-examined. 

Recent studies using modern radio line observations are 
consistent with no evolution in a^^„ii, a l beit w ith larger un- 
certainties than our result: I Wolfe etaH (1201 Ih obtain a 2a 
limit A{a^ g p fi) / g p jj. < 6 x 10"^ using 3C286 alone, based 
on the identification of the observ ed optical absorption with a 
specific H I-absorbing cloud, and lKanekar et alj (120101) com- 
pared 21 cm to UV C I absorption in two systems at z 1 .5 to 
find (6.8 ± l.Ostat ± 6.7sys) x 10"^. Recent UV-H I compar- 
isons obtaining similar-unc ertainty measurements to ours in- 
clude H^^^^^^^with (-0.1 ± 1.3) X lO"*" at (z) = 
1 .36 and lSrianandetaTI dloTol) who find (-1 .7 ± 1 .7) x lO"*" 
in a single absorber at z = 3.17. 



5. CONCLUSIONS 

We have obtained a precise constraint on the cosmic ac- 
celeration by directly measuring the real-time secular redshift 
drift of H 1 21 cm absorption line systems. We measure a red- 
shift drift of (z) = (-2.3 ± 0.8) X 10"** yr"' and an acceleration 



(v) = -5.5 ± 2.2 m s ' yr both consistent with no acceler- 
ation. These measurements have not yet reached the theoret- 
ically expected acceleration, z' = 2 x 10"^' yr"' or Av/Afo = 
0.3 cm s"' yr"' at z = 0.5, but they demonstrate an encouraging 
lack of peculiar acceleration in absorption line systems and 
the frequency stability of modern radio telescopes required 
for repeatable measurements. We thus find no significant sys- 
tematic effects that would impede even more precise measure- 
ments, which may be done by expanding the sample, using 
longer time baselines, larger redshifts, and higher signal-to- 
noise observations, particularly from larger apertures. Sys- 
tematic effects will be significant below 1 cm s"' yr"', but di- 
rect measurements of the cosmic acceleration using multiple 
methods, wavelengths, and telescopes will likely be possible 
in a few decades. 

We thank the GBT and NRAO staff for supporting this 
project and for building and operating an exquisitely sensi- 
tive and stable telescope. We also thank the referee for help- 
ful comments. This research made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, un- 
der contract with NASA. 



Amendola. L., Balbi, A., & Quercellini, C. 2008, Physics Letters B, 660, 81 
Balbi, A. & Quercellini, C. 2007, MNRAS, 382, 1623 
Briggs, F. H. 1983, ApJ, 274, 86 

Cai-illi, C. L., Menten, K. M., Stocke, J. T., Perlman, E., Vermeulen, R., 
Briggs, R, de Bmyn, A. G., Conway, J., & Moore, C. P. 2000, 
Phys. Rev. Lett., 85, 5511 

Corasaniti, P-S., Huterer, D., & Melchiorri, A. 2007, Phys Rev D, 75, 
062001 

Dai-ling, J., Macdonald, E. P, Haynes, M. P, & Giovanelli, R. 201 1, ApJ, 
742, 60 

Demleitner, M., Accomazzi, A., Eichhorn, G., Grant, C. S., Kurtz, M. J., & 
Munay, S. S. 2001, ASP Conf. Proc, Astronomical Data Analysis 
Software and Systems X, eds. F. R. Hamden, Jr, F. A. Primini, H. E. 
Payne, 238, 321 

Johnston, S., Taylor, R., Bailes, M., et al. 2008, Expeiimental Astronomy, 
22, 151 

Kanekar, N. & Chengalur, J. N. 2001, MNRAS, 325, 631 
Kanekar, N., Prochaska, J. X., Ellison, S. L., & Chengalur, J. N. 2010, ApJ, 
712, L148 

Lane, W. M. & Briggs, E H. 2001, ApJ, 561, L27 

Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 

Liske, J., Grazian, A., Vanzella, E., et al. 2008, MNRAS, 386, 1 192 



Loeb, A. 1998, ApJ, 499, LI 11 

Rahmani, H., Srianand, R., Gupta, N., Petitjean, P., Noterdaeme, P., & 

Vasquez, D. A. 2012, MNRAS, 425, 556 
Sandage, A. 1962, ApJ, 136, 319 

Srianand, R., Gupta, N., Petitjean, P., Noterdaeme, P., & Ledoux, C. 2010, 

MNRAS, 405, 1888 
Titov, O., Lambert, S. B., & Gontier, A.-M. 2011, A&A, 529, A91 
Tzanavaris, P, Murphy, M. T., Webb, J. K., Flambaum, V. V., & Curran, S. J. 

2007, MNRAS, 374, 634 
Webb, J. K., King, J. A., Murphy, M. T., Flambaum, V. V., Carswell, R. F, & 

Bainbridge, M. B. 2011, Phys. Rev. Lett., 107, 191101 
Wolfe, A. M., Brown, R. L., & Roberts, M. S. 1976, Phys. Rev. Lett., 37, 179 
Wolfe, A. M., Broderick, J. J., Johnston, K. J., & Condon, J. J. 1978, ApJ, 

222, 752 

Wolfe, A. M., Davis, M. M., & Briggs, F H. 1982, ApJ, 259, 495 

Wolfe, A. M., Jorgenson, R. A., Robishaw, T., Heiles, C, & Prochaska, J. X. 

2011, ApL 733, 24 
Xu, M. H., Wang, G. L., & Zhao, M. 2012, A&A, 544, A135 
Zakamska, N. L. & Tremaine, S. 2005, AJ, 130, 1939 



